Proceedings of the DPF-2009 Conference, Detroit, MI, July 27-31, 2009 



1 



Model-Independent Results for the Decay B+ i^v'i'y at H4B>V7 

D. M. Lindemann (On behalf of the BaB^r Collaboration) 

Department of Physics, McGiil University, Montreal, Quebec, Canada H3A 2T8 



O 

o 

(N 
Oh' 

O 

m 



X 
CD 

Oh 

CD 



> 

ON 
00 

\o 
vn 

a\ 
o 

ON 

o 



X 



We present a search for the radiative leptonic decays _B+ — ► e+i>e7 and -B"*" — > using data collected by 

the BAB^, detector at the PEP-II B factory. We fully reconstruct the hadronic decay of one of the B mesons 
in T(45) — > B^B~ and then search for evidence of the signal decay within the rest of the event. This method 
provides clean kinematic information on the signal's missing energy and high momentum photon and lepton, 
and allows for a model-independent analysis of this decay. Using a data sample of 465 million B-meson pairs, 
we obtain sensitivity to branching fractions of the same order as predicted by the Standard Model. We report 
a model-independent branching fraction upper limit of B{B+ i+vri) < 15.6 x 10"'' (£ = e or ^i) at the 90% 
confidence level. 



1. Introduction 

The leptonic decay i'^i^g'y, where £ = e or fi, 

proceeds via an annihilation of b and u quarks into 
a virtual boson with the radiation of a photon.^ 
Leptonic decays can provide clean theoretical predic- 
tions of Standard Model (SM) parameters without the 
QCD-based uncertainties arising from hadrons in the 
final state. The purely leptonic decay i^i'e, 
which offers a clean prediction of the _B-meson decay 
constant fs, is helicity suppressed, having a branch- 
ing fraction that is proportional to the square of the 
lepton mass. Although the radiative mode is ad- 
ditionally suppressed by a factor aom, the presence 
of the photon can remove the helicity suppression 
by affecting the coupling of the spin-0 B meson to 
the spin-1 boson, possibly through an interme- 
diate off-shell state [l|. The branching fraction of 
5+ i^i'e'y is predicted in the SM to be indepen- 
dent of the lepton type and of order 10~^, making it 
potentially accessible at current and future B facto- 
ries such as BaBAR and SuperB ^ respectively. The 
most stringent published limits are from the CLEO 
coUaboration with B{B^ e+i/gl) < 2.0 x 10"'' and 
B{B+ < 5.2 X 10-^ at the 90% confidence 

level (CL) ^. 

The branching fraction of the signal decay can be 
written as Q: 



B{B+^ i+vn) = 



2887r^ 



I Kb I fB^^TB \- 



A, 



(1) 

where Gi? is the Fermi constant, Vub is the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix element describ- 
ing the coupling of b and u quarks, friB and tb are 
the S-meson mass and lifetime respectively, Qi is the 
quark charge, and is the first inverse moment of 
the light-cone i3-meson wave function. This last pa- 



^ Charge conjugate modes are included implicitly throughout 
this paper. 



rameter plays an important role in proving QCD fac- 
torization [5|. It also enters into calculations of the 
i? — > TT form factor at zero momentum transfer and the 
branching fractions of two-body hadronic i3-meson de- 
cays such as B ^ Dtt and B ^ nn \d\, the latter be- 
ing a benchmark channel for measuring the angle a 
of the CKM Unitarity Triangle. Typically assumed to 
be of order Aqcd, at a few hundred MeV, Xb cur- 
rently suffers from significant theoretical uncertainty 
0. A measurement of B{B^ — > i^i^el) can provide a 
clean prediction of this important parameter. In ad- 
dition, since B^ i^vij is a possible background for 
^ a measurement of B{B^ — I'^vi'j) over 

the full photon energy spectrum is needed for an accu- 
rate measurement of B{B^ I'^vi) [§], which would 
improve predictions of /s. 



2. Analysis 

This analysis uses data taken from the BaBAR exper- 
iment 0, located at the asymmetric-energy PEP-II 
e+e" storage rings at SLAC in California. On-peak 
data is produced by colliding electrons and positrons 
with a center-of-mass (CM) energy at the T{AS) reso- 
nance of 10.58 GeV, just above the threshold for BB 
production. We use the full BABMt dataset of 465 ± 5 
million i?-meson pairs, corresponding to an integrated 
luminosity of 423 fb^^. Charged-particle momenta 
are measured using a tracking system comprised of 
a silicon vertex detector and drift chamber contained 
within the uniform magnetic field of a 1.5T super- 
conducting solenoid. Charged-particle identification 
is based on the energy loss in the tracking system 
and the Cherenkov angle in a ring-imaging Cherenkov 
detector. Photon energies and electron identification 
are provided by a CsI(Tl) scintillating electromagnetic 
calorimeter (EMC). Finally, muons are distinguished 
from hadrons using instrumentation embedded within 
the steel magnetic-flux return. Monte Carlo (MC) 
simulations are used to model the detector response, 
determine the signal efficiency, and study the back- 
ground. 
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2.1. Hadronic Reconstruction 

We present the first search for B+ t'^vg-y that 
uses an exclusive method, by fuUy reconstructing one 
B meson (Stag) via hadronic decay modes. Unlike 
previous inclusive searches of B+ — > t^vg'^, which 
identified the signal lepton and photon candidates be- 
fore the recoiling B meson, we reconstruct the re- 
coil B meson first before searching for the signal de- 
cay. Although this technique results in a low sig- 
nal efficiency (0.3% for signal modes), and thus more 
statistically-limited results, it compensates by provid- 
ing a highly pure sample of B mesons with compar- 
atively little non-BB (continuum) background. Our 
reliance on poorly-modeled continuum background is 
thus reduced, which enables a search for i?+ — *■ t^vti 
that can avoid kinematic constraints in the signal 
selection. Therefore, we present the first model- 
independent search of B(B^ i^vi'j) over the full 
kinematic range. In addition, by reconstructing the 
Stag using only detectable hadronic decay modes, the 
missing four- vector of the otherwise undetectable sig- 
nal neutrino is fully determined. 

The event selection begins by reconstructing a Z)'*-' 
meson using various hadronic decay modes, listed in 
Ref. . Additional charged tracks and neutral EMC 
clusters are then combined with the D'^*^ seed to form 
the decay B —> _D'^*^Xhad, where Xhad is a combina- 
tion of neutral and charged pions and/or kaons. Xhad 
is chosen such that \/\E\ = l-Bst.g - ^1 < 0.12, where 
^/s is the total energy of the e+e" system and Sst^g is 
the Stag candidate energy, both in the CM frame. If 
there exists more than one Stag candidate, we choose 
the one with the D^*^ and Xhad mode that provides 
the highest purity of well-reconstructed candidates. 
Finally, we require that the Stag candidate is charged. 

Since the beam energy at PEP-II is accurately 
known within a few MeV, we define the mass of the 

Stag as mEs = ^1 "P-Btag' where pb^^^ is the Stag 
three-momentum in the CM frame. Since the Stag 
candidate should have a mass consistent with the 
nominal S-mcson mass, we require 5.27 < toes < 
5.29GeV/c^. The toes sideband region, defined as 
5.20 < Toes < 5.26GeV/c^, is populated by combina- 
toric background, events coming from a T{AS) BB 
or continuum decay in which the Stag is incorrectly 
reconstructed. To improve the MC estimation of the 
reconstruction efficiency, we normalize the MC that 
peaks within the toes signal region to the data that 
peaks within this same region, after applying the Stag 
selection criteria. In addition, we reduce our reliance 
on the MC estimate of combinatoric events by esti- 
mating the non-peaking background directly from the 
data within the toes sideband. 

When a T(4S') resonance decays to a BB pair, the 
two particles are almost at rest, with a momenta of 
about 350MeV/c in the CM frame, since their masses 



are about half the T{4S) mass. Therefore, they tend 
to decay with an isotropically symmetric topology in 
the CM frame. On the other hand, lighter qq and 
r+r~ pairs, which are also produced in the e^e~ 
interactions, have a higher momentum in the CM 
frame. Their decays tend to have a more jet-like 
shape with a strongly-preferred direction characteriz- 
ing the event, preferentially at small angles in relation 
to the beam axis. Therefore, to achieve a higher pu- 
rity of S mesons, we use a multivariate selector of five 
event-shape variables to separate BB from continuum 
events. We require: 



> 30%, (2) 



where Vsixi) and Vq{xi) are probability density func- 
tions determined from MC that describe BB and con- 
tinuum events, respectively, for the five event-shape 
variables x;. These variables, as explained in Ref. [lO|, 
describe the sphericity of the event, the magnitude 
and directions of the thrust axes, and the direction of 
the Stag momentum with respect to the beam axis. 
This requirement also suppresses the poorly-modeled 
continuum backgrounds which can contribute to dis- 
crepancies between the MC and data. 



2.2. Signal Selection 

After the Stag reconstruction, we assign all the re- 
maining charged tracks and neutral EMC clusters, to 
the signal S meson (Sgig), as well as any missing mo- 
mentum (piniss) within the event. We require that Sgig 
has exactly one track, with a charge opposite that of 
the Stag charge. In addition, this signal track must 
satisfy particle identification criteria of either an elec- 
tron or muon and fail that of a kaon. 

Because high energy electrons can emit 
bremsstrahlung photons, we search for signal- 
side clusters that have proximity to the EMC deposit 
of an electron candidate, both in 6, the angle from 
the beam axis, and in (p, the azimuthal angle around 
the beam axis. Because the detector has a solenoid 
magnet through which all the particles travel, charged 
particles are bent in (f) according to their charge, and 
therefore is multiplied by the lepton candidate 
charge. Any cluster with a momentum vector that 
is separated from that of the electron-identified 
signal track by |A6'| < 3° and -3° < < 13° is 
identified as a bremsstrahlung photon. The energy 
of this cluster is then used to correct the four-vector 
of the signal track [pi) and is removed from the 
list of clusters. Next, we identify the signal photon 
candidate, whose energy spectrum is expected to 
peak around 1 GeV. We do this by searching through 
the remaining signal-side clusters for the one with 
the highest energy in the CM frame. 
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Although 5+ t^vg'j should only have one signal- 
side cluster if correctly reconstructed, extra clusters 
within an event can be due to fragments from particle 
showers, low-energy clusters from the -Stag, and/or 
noise in the detector such as beam-related photons. 
Because of the presence of extra clusters, we apply the 
very loose requirement that there are no more than 11 
additional low energy clusters which are not otherwise 
used in the i?tag or signal reconstruction. The total 
energy of these clusters is required to be less than 
800 MeV, where we only include in this sum clusters 
which have lab-frame energy greater than 50 MeV. 

Missing energy can be due to either an undetectable 
particle such as a neutrino or a detectable particle 
that travels outside of the fiducial acceptance of the 
detector. To suppress events of the latter category, we 
require that Pmiss points within the fiducial acceptance 
of the detector. 

The lepton and neutrino are produced back-to-back 
by the virtual boson. However, in the Bsig rest 
frame, the angle between the signal track momen- 
tum and Pmiss is affected by the release of the photon. 
Therefore, we require cosOii, < —0.93 in the frame re- 
coiling from the photon. This frame is defined as the 
difference between the B^ig four- vector (ps) and the 
signal photon candidate four- vector (p-^). Including 
this requirement assures that the decay is consistent 
with being a three-body decay, independent of when 
the photon was released. 

The most discriminating variable in this analysis is 
a requirement on the neutrino candidate's invariant 
mass m^, given by: 



{PB - P7 - Pi) 



(3) 



This requirement isolates events in which the pho- 
ton and lepton candidates are in kinematic agree- 
ment with a massless third daughter in the event. 
Fig. [T] shows that the signal peaks at zero, while the 
the background rises with m^. The tail on the sig- 
nal distribution is larger for the electron mode than 
the muon mode due to unrecovered bremsstrahlung 
photons. The signal region for this variable is de- 
fined as -1 < m2 < 0.46 (0.46) GeVVc'* for the elec- 
tron (muon) mode. These values, as with all other cut 
values in this analysis, were optimized using the figure 

) [11], where 



of merit £^^^ /{\ 



N, 



bkg 



bkg 



1.3, e) 



is the number 



is the total signal efficiency, and iV^ 
of expected background events. 

The and cos^^jy requirements kinematically re- 
strict the types of events that pass. Since light 
mesons often decay to a pair of photons, where Xu is a 
neutral meson containing a w-quark, the topology and 
kinematics of a B'^ — > X^^l^vi event can often mimic 
that of the signal decay. This is particularly true if one 
photon is "missing" such as through misreconstruc- 
tion into the Stag or due to a lab-frame energy that is 
too low to be detectable by the EMC. Thus, the pri- 



> 

CM 

d 



> 
d 




'2^(GeV2/c4) 



Figure 1: distribution after all selection criteria are 
applied, in electron (top) and muon (bottom) modes for 
the mES-peaking (shaded) plus non-peaking (solid) con- 
tributions in the full background MC sample, signal MC 
normalized to B = 40 x 10~® (dashed), and data (points). 
Events to the left of the vertical lines are selected. 



mary background that passes the signal selection are 
_B+ — > X^l^vi events with a high energy photon. To 
suppress B^ X^f+i/^ events in which both photon 
daughters are present in the signal-side clusters, we re- 
ject events in which the signal photon candidate can 
be combined with another cluster to form an invariant 
mass consistent with the tt" or rj mass. To improve 
the purity of the reconstructed tt*' mesons, we only 
use clusters that are above a given energy Ej2 in the 
-Bsig rest frame. To suppress B~^ 'n'^t^vn events, we 
reject events in which the invariant mass is between 
120 and 145 MeV with > 30 MeV or between 100 
and 160 MeV with i?^2 > 80 MeV. Since 77 particles are 
much more massive and thus tend to decay into higher 
energy photons, we suppress B^ rjl^v^ events by 
rejecting events with E^2 > lOOMeT^ and an invariant 
mass between 515 and 570 MeV. Finally, we suppress 
B^ uj£~^i^g [tt'^ j]£~^ vi events by rejecting events 
in which the signal photon candidate combined with a 
TT*^ candidate forms an invariant mass between 115 and 
145 MeV. This 7r° candidate must have E^2 > 70 MeV 
and an invariant mass between 115 and 145 MeV. 

Events in which the two photons from a, B^ 
7r°£+i/f decay are merged into a single EMC cluster 
can mimic the signal kinematics exactly, since the 
cluster chosen as the signal photon candidate contains 
the full energy of the 7r°. We suppress this background 
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using a requirement on the cluster width of the signal 
photon candidate. A cluster that is the result of two 
merged photons from a tt" decay tends to be wider 
than a cluster from a single photon. We require the 
lateral moment to be less than 55%. 



2.3. Backgrounds and Uncertainties 

We define the signal branching fraction for each lep- 
ton mode t as: 



bkg 



(4) 



where is the number of observed data events 

within the signal region and Ng± = 465 x 10^ is 
the number of mesons in the data sample. The 
branching fractions are computed using the frequen- 
tist formalism of Feldman and Cousins [l^ , with un- 
certainties on N^^^ and e^'^ modeled using Gaussian 
distributions. Since B{B^ t'^vg'j) is expected to be 
independent of the lepton type, we also combine the 
branching fractions of the two modes. 

The MC indicates that the selection criteria, es- 
pecially the one track requirement and the and 
cos 0^1, restrictions on the event kinematics, removes 
all decay types except X^l^vg events and 

the occasional combinatoric event from a misrecon- 
structed -Stag- Thus, we split N^^^ into two cat- 
egories: peaking events {Nj°^^) and combinatoric 
events {Nj""^""). An event is considered peaking if 
it comes from a T{AS) — > BB decay in which the 
-Btag is correctly reconstructed and hence peaks within 
the TTiEs signal region. Since only i?+ X^t^Vf^ 
events peak in this region, we improve the statistics of 
by using exclusive -B+ ^ X^l'^vi MC samples. 
For Nf°™^ , we extrapolate the number of misrecon- 
structed -Btag events in the mss signal region directly 
from the number of data events in the tties sideband. 

The uncertainty on N'^°'"^^ is dominated by the side- 
band data statistics, but it also includes a 14.6% sys- 
tematic uncertainty on the toes combinatoric back- 
ground shape. On the other hand, since e^'® and 
Nf^^ rely on the accuracy of the MC simulations, 
their uncertainties arise not only from MC statistics, 
but also include systematic uncertainties due to how 
well the MC agrees with the data. The contributions 
to the systematic uncertainties are listed in Table HI 
The systematic uncertainty of Nf'^'^^ is dominated by 
a 13.6% uncertainty in the branching fractions and 
form factors of various exclusive 5+ X'^£~^V£ decays 
[Tot . The uncertainty due to the Stag reconstruction, 
which also accounts for the uncertainty in , is de- 
termined by varying the shape and scaling of the toes 
combinatoric distribution. 



Table I Contributions (in percent) to the systematic un- 
certainty on the branching fraction due to the signal ef- 
ficiency e"^ and the peaking-background estimate Nf"^^. 
The B{B+ X°e+ve) source refers to branching fraction 
and form factor uncertainties in — > XHi'^ug decays. 



Source of 
systematics 


e mode 


ig 

/i mode 


e mode 


Gak 

/i mode 


B{B+ X^i+i^i) 






13.6 


13.6 


-Btag reconstruction 


3.1 


3.1 


3.1 


3.1 


Particle identification 


0.9 


1.3 


0.9 


1.3 


Track reconstruction 


0.4 


0.4 


0.4 


0.4 


Photon reconstruction 


1.8 


1.8 


1.8 


1.8 


Cb 


1.4 


1.4 


1.4 


1.4 


ml 


0.5 


0.5 


1.4 


1.4 



3. Results 

To avoid experimenter bias, we blinded the data in 
the region ml < 1 GeV^/c^ and toes > 5.26GeV/c^ 
until all selection criteria, background estimates, and 
systematic uncertainties were finalized. The final sig- 
nal efficiencies correspond to approximately one signal 
event per mode, assuming a branching fraction within 
the SM predictions. After unblinding, we observe 4 (7) 
data events within the signal region for the electron 
(muon) mode, compared to an expected background 
of 2.7 ± 0.6 (3.4 ± 0.9) events. The branching fraction 
results are given in Table HIl 



rbkg 



^comb ^ 



Table 11 Expected background yields N^. = iti 
Nf'^^'^ , signal efficiencies e^'^, number of observed data 
events N^^", resulting branching fraction limits at 90% 
CL, and the combined central value Bcomtined- Uncertain- 
ties are given as statistical ± systematic. 





B+ e+i^e7 


B+ ^i+i^^-f 


jycomb 
,rpcak 


0.3 ± 0.3 ± 0.1 
2.4 ± 0.3 ± 0.4 


1.2 ± 0.6 ± 0.6 
2.1 ± 0.3 ± 0.3 


Big 
'=■1 


2.7 ± 0.4 ± 0.4 
(7.8 ± 0.1 ± 0.3) X 10"* 
4 


3.4 ± 0.7 ± 0.7 
(8.1 ± 0.1 ± 0.3) xlO^* 
7 


B combined 


(6.5lI;?l-«) X 10-« 


Model-ind. 
Limits 


<17xl0"'^ 1 < 26x10"'' 
< 15.6 X 10"*' 



Although the effective detector and particle identi- 
fication thresholds are about 20MeV for photon en- 
ergy and 400 (800) MeV for electron (muon) momen- 
tum, and we apply no minimum energy requirements. 
Thus, this analysis is essentially valid over the full 
kinematic range, as shown in Fig. [2l However, there 
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exists theoretical uncertainty in the photon energy 
spectrum below Aqcd for the — > i^vgj decay. 
Therefore, using certain theoretical techniques, the 
extraction of As can be improved by including a min- 
imum energy requirement on the signal photon ^TJ. 
When the signal photon candidate energy is required 
to be greater than 1 GeV, we observe 2 (4) data events 
with N^^'^ = 1.4 ± 0.3 (2.5 ± 1.0) in the electron 
(muon) mode. The e^'^, which is mostly uncorre- 
lated with the photon energy spectrum, is reduced 
by 30%, resulting in a partial branching fraction of 
AB{B+ i+va) < 14 X 10-6 at 90% CL. 

The differential branching fraction versus photon 
energy of i?"*" — > i'^vci is given by: 



487r2 



(5) 

where y = 2E^/mB- The two form factors, /y and /a, 
describe the vector and axial- vector contributions, re- 
spectively, to the i3 ^ 7 transition. Although /a = 
in some models ^Ij and was assumed in the CLEO 
measurement of B{B^ i^i^ij), most models assert 
/a = /y [ll- In this analysis, we use signal MC that is 
generated based on the tree-level hadronic matrix el- 
ement for 5+ €^f£7 as described by Ref. Q, using 
a minimum photon energy of 350 MeV. We use both 
form-factor models for our signal MC to evaluate the 
impact of the decay model on the signal selection effi- 
ciency and to ensure model-independency. We deter- 
mine e^g^ using the /a — fv signal model, but because 
our analysis is independent of the decay kinematics, 
the /a = model yields consistent ef^ values. 

We also determine branching fraction limits that are 
dependent on the signal model by introducing a kine- 
matic requirement on the angles between the three 
daughter particles of the signal decay. We use the 
quantities cos^-y^ and cosO^jj, where O^i is the an- 
gle between the photon candidate and signal track 
momenta, and 0^^, is the angle between the photon 
candidate momentum and Pmiss, both in the Bgig rest 
frame. As seen in Fig. [31 the photon is emitted prefer- 
entially back-to-back with the Icpton in the /a — fv 
model, and back-to-back with either the lepton or 
neutrino in the /a = model. Thus, we require 
(cos 6*^^.- 1)2 + (cos 0^,. -I- 1)2/3 > 0.4 or (cos6'y^-l)2-h 

: model, and only 
fv- This reduces 
e"^ in both modes and models by 40%. We observe 
(0) data events in the electron (muon) mode with 
N^^^ = 0.6 ± 0.1 (1.0 ± 0.4) for the /a = fv model, 
corresponding to B{B^ (.^vi'^) < 3.0 x 10^^. Like- 
wise, in the /a = model, we observe 3 (2) data 



(008 6*^^ -I- 1)2/3 > 0.4 for the /a 
the former relationship for /a 



events with N^^^ = 1.2 ± 0.4 (1.5 ± 0.6), correspond- 
ing to B{B+ i+va) < 18 X 10-6. 

4. Conclusion 

In conclusion, we have searched for i?+ —t t^vci 
using a hadronic recoil technique and observe no sig- 
nificant signal within a data sample of 465 million 
BB pairs. We report a model-independent limit of 
B{B^ -> t^va) < 15.6 X 10-6 at the 90% CL, which 
is consistent with the standard model prediction and 
is the most stringent published upper limit to date. 
Using Eq. (P) with fs = 0.216 ± 0.022 GeV, ms = 
5.279 GeV/c2, tb = 1.638 ps, mb = 4.20GeV/c2, and 
\Vub\ = (3.93±0.36) X 10-3 [l3|, the combined branch- 
ing fraction likelihood function corresponds to a limit 
of As > 0.3 GeV at the 90% CL. 
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Figure 2: The distribution of the signal track momentum versus the signal photon candidate energy, in the Bsig rest 
frame, after all signal selection criteria is applied, for data (points), /a = fv signal model (left), /a — O signal model 
(middle), and N^^^ (right). The electron modes are shown on the top and the muon modes are on the bottom. The 



size of the box in the A'^^''^ plots is proportional to the number of background events within the histogram bin. The 



contribution of Np 



which is determined using data events in the tties sideband, is shown as red boxes. 




Figure 3: The correlation between cos^^.^ and cos6^„, in the _Bsig rest frame, after all signal selection criteria is applied, 
for data (points), /a = fv signal model (left), fA=0 signal model (middle), and A'^''^ (right). The electron modes 
are shown on the top and the muon modes are on the bottom. The size of the box in the N^^^ plots is proportional 
to the number of background events within the histogram bin. The contribution of N^"'^^, which is determined using 
data events in the tties sideband, is shown as red boxes. The black arcs indicate the cut-off values for the model-specific 
requirements. 



